Observation of Giant Positive Magnetoresistance in a Cooper Pair Insulator 
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Ultrathin amorphous Bi films, patterned with a nano-honeycomb array of holes, can exhibit an insulating 
phase with transport dominated by the incoherent motion of Cooper Pairs (CP) of electrons between localized 
states. Here we show that the magnetoresistance (MR) of this Cooper Pair Insulator (CPI) phase is positive and 
grows exponentially with decreasing temperature, T, for T well below the pair formation temperature. It peaks 
at a field estimated to be sufficient to break the pairs and then decreases monotonically into a regime in which 
the film resistance assumes the T dependence appropriate for weakly localized single electron transport. We 
discuss how these results support proposals that the large MR peaks in other unpatterned, ultrathin film systems 
disclose a CPI phase and provide new insight into the CP localization. 



Below its transition temperature, T c0 , a conventional super- 
conductor, like Pb, can be driven into its non- superconducting, 
normal state by applying a magnetic field, H. The tempera- 
ture, T, dependent sheet resistance, R\j(T), of this state joins 
smoothly to the normal state resistance just above T c o, Rn, 
and assumes the dependence expected for a simple metal of 
unpaired electrons Q. The behavior of this low T normal state 
changes substantially when these superconductors are made 
as thin films with Rn ~ Rq = h/Ae 2 . For example, apply- 
ing a H to superconducting (SC) films of either Indium oxide 
(InO^ ) O SI or TiNl5l drives them directly into an insulat- 
ing phase. Their R\j (T) rise exponentially with decreasing T 
to exceed .Rat by orders of magnitude. A similarly dramatic 
Superconductor to Insulator Transition (SIT) can be achieved 
in granular Pb films by increasing RnM&- These behaviors 
have been taken to imply that the films enter a Cooper Pair In- 
sulator (CPI) phase consisting of exponentially localized, but 
intact Cooper Pairs of electrons (3) O [TJ [8) . 

For granular Pb films, the formation of a CPI phase has 
strong intuitive appeal and experimental support. STM exper- 
iments show that they consist of islands of grains that can nat- 
urally localize CPs|9|. Indeed, tunneling experiments on in- 
sulating films confirmed the existence of these localized pairs 
by showing the energy gap in the density of states that accom- 
panies CP formation! 6]. Also, these insulators exhibit giant 
negative Magneto-Resistance (MR) that can be attributed to 
the enhancement of inter-island quasi-particle tunneling! 10]. 
By contrast, InO^ and TiN films lack any obvious structure 
that could localize CPs. Moreover, these films exhibit a gi- 
ant positive MR0 |U [TTJ [El, which can peak orders of 
magnitude above Rn at sufficiently low T. The mechanism 
behind this spectacular giant positive MR0 ffl |5j [TTJ [121 and 
whether it is a property of a CPI phase remains unresolved 
despite significant attention Hl3\ [T4l [T5l [T6l . 

Theories of the positive MR peak presume that CPs spon- 
taneously localize into islands or puddles fT3l 1T41 [T5l IT6l . On 
each island, the complex SC order parameter has a well de- 
fined amplitude, but electrostatic interactions between islands 
prevents the development of the long range phase coherence 
necessary for CP derealization. A magnetic field induces 
more phase disorder and localization through the direct cou- 




Figure 1 : a) SEM image of the nano-honeycomb substrate. The hole 
center to center spacing and radii are 100 ± 5 and 27 ± 3 nm, re- 
spectively. Arrows denote H. b) Sheet resistance as a function of 
temperature, i?n (T), of NHC films produced through a series of Bi 
evaporations. The film 16 is the last insulating film and SI is the first 
superconducting film in the series, c) Surface plot of i?n (T, H) for 
film 16, which has a normal state sheet resistance of 19.6 kQ and 1.1 
nm Bi thickness. The solid lines are isotherms. 



pling of the vector potential to the order parameter phase and 
by reducing the order parameter amplitude through its CP 
breaking effects. At very high H, the latter effect leads to 
a negative MR as the film returns to an unpaired state lfTTl . 
Surprisingly, this basic picture qualitatively accounts for the 
positive MR peak of smooth InO^ and TiN films. This agree- 
ment and recent STM results fT8lL suggest that a CPI phase 
can spontaneously develop due to electronic interactions |[T9lL 
H\20\ and/or disorder induced localization [21 j. Verifying the 
connection between the MR peak and the CPI phase requires 
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creating a system exhibiting similar MR features that can be 
simultaneously probed for the CPI phase. 

Here we present the MR of an amorphous film system with 
a clear CPI phase. These ultrathin Bi films are patterned with 
a Nano-HoneyComb (NHC) array of holes[7|. They exhibit a 
large positive MR that peaks at a field comparable to the es- 
timated average depairing field. Superimposed on the initial 
rise of the resistance are oscillations with a period set by the 
SC flux quantum indicating that localized CPs dominate the 
transport. Moreover, the R\j (T) are activated and the oscil- 
lations and the rise in R\j (H) stem primarily from variations 
in the activation energy. At fields well beyond the peak, the 
transport appears dominated by weakly localized quasiparti- 
cles. Many of these characteristics of the MR of NHC films 
are shared by unpatterned InO^ and TiN films, which sup- 
ports earlier proposals that their peaks reflect an underlying 
CPI phase with transport dominated by CP motion. Further- 
more, the behavior of the activation energy of NHC films is 
consistent with H controlling the CP localization by tuning 
the energy characterizing the coupling of the relative phases 
of the localized CP states. 

The data presented characterize amorphous, NHC Bi films 
that were deposited in situ and measured in a custom designed 
dilution refrigerator cryostat|8|. The films were patterned by 
depositing them onto aluminum oxide substrates structured 
with a NHC array of holes (see Fig. [TJ[ 22 1. Repeated Bi 
evaporations yielded a series of films spanning the thickness 
tuned SIT (see Fig[IJ))|7'|. Film sheet resistances as a function 
of T and H, R\j(T, H), were measured using standard four- 
point AC and DC techniques employing currents sufficiently 
small to be in the linear regime of the current- voltage charac- 
teristics. A SC solenoid provided H up to 8T transverse to the 
film plane. The data presented here come from the two films 
nearest the SIT, 16 and S 1 . Experiments on films on two other 
substrates yielded very similar results. 

Patterning these Bi/Sb films makes it possible to detect the 
charge of the carriers involved in their transport!?) and ap- 
pears to be essential to the formation of the CPI phasef7ll23l. 
The hole array imposes a spatial period with a unit cell area S, 
which sets the H scale, Hm = ®o/S, corresponding to one 
SC flux quantum, $0 = h/2e, per unit cell of the array. As 
with Josephson Junction and wire arrays, the thermodynamic 
and transport properties of NHC films are expected to oscillate 
with period HmTlI. We will use the frustration, / = H/Hm 
to specify H with H M = 0.21T. 

The MR of NHC films near their thickness-tuned SIT ex- 
hibit a rich structure including oscillations at low H that 
merge into a giant peak at higher H. The i?n(T, H) surface 
shown in Fig. [T] demonstrates these behaviors for the insu- 
lating film, 16 (see Fig. [lj>). The surface was generated by 
interpolating data acquired by sweeping T at constant H for 
a series of closely spaced H. At low H and T, the MR oscil- 
lates with the period expected for CPs indicating that the film 
is in the CPI state|7|. The oscillations first appear for T ~ 
0.6K, which gives a lower limit for the CP formation temper- 
ature, T c q. At the lowest T, five oscillations are easily re- 
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Figure 2: Sheet resistance as a function of if at 100 and 120 mK for 
films 16 and SI. The lines are spline fits to the data points (shown as 
symbols on the 100 mK traces). Inset: Magnified view of the low H 
data. 



solved and there is a hint of a sixth. Oscillations 4-6 are more 
difficult to resolve because they diminish in amplitude and ap- 
pear on a rapidly rising background MR. At 100 mK, R\j (H) 
peaks near f = 7 with a resistance that is nearly a factor of 
30 larger than the zero field value. This giant MR peak grows 
and moves to lower field with decreasing T. These R(T, H) 
features evolve smoothly with decreasing Rn and across the 
thickness tuned SIT as shown by Fig. [2] . The peak moves to 
higher H and diminishes in size and the oscillation amplitude, 
measured on a linear scale, shrinks. We hasten to add that the 
substrate patterning! 23 ] is essential to the appearance of the 
MR peak. Similar Bi/Sb films deposited on glass substrates 
show only a very small (< 10%) peak[24 |. 

In the low T limit, the R\j (T) fit well to an Arrhenius form, 
Rn(T) = Roe T °/ T , for H values extending to the peak (see 
Fig. [3ji). While only R n (T) at integer / are shown, R n (T) 
at non-integer / behave similarly. The activation energy, To 
qualitatively mirrors the MR showing large oscillations with 
the same period and climbing to a maximum near / = 6 of 
To = 0.35 K (see Fig. Moreover, the peak To decreases 
to 0.27 K and moves out to / ~ 9 for the lower R^ film, SI 
(not shown). The prefactor, Ro, grows relatively slowly up to 
the peak field (see Fig. [3}:). Its factor of 2 increase is substan- 
tially smaller than the size of the MR peak at 100 mK. Taken 
together, these behaviors indicate that the MR oscillations and 
the peak in the MR stem primarily from an H dependence of 
the activation energy. 

Well beyond the peak (jjlqH > 4T), where the R(H) 
appear to asymptotically approach a T dependent constant, 
the conductance, G(T) = iT^T) fits better to G(T) = 
AG oln(T/T ) + G(T ) than an activated form. Fig. [^com- 
pares these forms for 16 and SI at 6.8 T. The activated form 
can only account for the data over a limited range, while the 
logarithmic dependence extends over nearly a decade. This 
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Figure 3: a) The solid lines are the sheet resistance versus inverse 
temperature for film 16 in fields corresponding to / = 0,1,2,3,4,5,6,7. 
The dashed lines are fits to R n (T) = R exp(T /T). b) and c) 
exhibit the fitting parameters To and Ro, respectively. 



4r 

3 

P 

1*2 



r 



r 



0.5 1 

T C0 (K) 



1.5 



Figure 5: Magnetic field at which the activation energy peaks versus 
the temperature at which R(T) has dropped to 0.8 of its maximum. 
The data points were acquired from films deposited on three different 
substrates. Similarly shaped points (open and closed) come from 
adjacent squares of the same film. Films 16 and SI are the 6-pointed 
star and circle, respectively. The dashed line is a weighted linear fit 
to these data, whose slope is 1.7±0.2 T/K and y-intercept=0.3 T. 
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Figure 4: Film conductance, G(T) = 1/Rq normalized by the 
quantum of conductance, Goo = e 2 /27v 2 h versus temperature on a 
logarithmic scale for films 16 and SI in H = 6.8 T. Inset: Arrhenius 
plots of the same data shown in the main plot. 



fit, with A = 0.8 ~ 1, is consistent with transport by unpaired 
electrons in the presence of disorder enhanced e-e interactions 
in two dimensions [25]. 

The oscillations on the low field side of the MR peak inti- 
mate that CP transport processes are at the heart of the giant 
MR. For film 16, the oscillations extend to within 20% of the 
peak resistance. Interestingly, SI displays a similar number 
of oscillations as 16, but the oscillations do not appear as high 
on its peak, which appears at a greater field. This comparison 
suggests that the number of oscillations is limited by factors 
such as disorder in the hole array |'26'| rather than the physical 
effects giving rise to the peak. Thus, it is possible that CPs 



exist up to and beyond the peak. 

The values of the peak activation energy, T™ ax ', and mag- 
netic field, Hp, support the picture of CP dominated transport 
at low H and the emergence of single electron transport at 
high H. First, Tq 1 ^ consistently falls below the CP binding 
energy, 3.5/c£T c o, indicating that the low H transport does not 
involve the breaking of CPs] 17]. Second, H p appears related 
to the upper critical H, H C 2, the field required to break CPslfTTl. 
For amorphous films of constant resistivity, dirty limit theory 
implies, H C 2 cx T c ol 1 ]. We estimate T c o for the NHC films as 
the T at which R(T) has dropped to 0.8 of its maximum. We 
attribute this drop to the onset of very strong pairing fluctua- 
tions and fluctuation paraconductivity] 1 ]. The choice of 0.8 
rather than the standard 0.5 allows analysis of both SC and 
reentrant films. Other choices from 0.75 to 0.95 yield simi- 
lar results for the slope, but cause the intercept to vary. Data 
from six films from three different substrates with resistivities 
estimated to be similar to about 20% indicate that H p rises ap- 
proximately linearly with T c0 (see Fig|5]) with a slope of 1.7 
±.2 T/K. This slope is somewhat larger than what can be es- 
timated from previous experiments on unpatterned films (23, 
which implied that H c2 = 1.2T c0 . Thus, T™ ax appears at or 
just above the H necessary to depair most of the electrons. 

These giant MR peaks closely resemble those seen in un- 
patterned InO^H CED CE2 and TiNE [[6) films. Those also 
grow exponentially with decreasing T(31 [TT] [T6), become 
stronger with disorder, occur near the estimated H C 2\ 11 ], and 
evolve smoothly across the disorder tuned SIT[5|. In addi- 
tion, their transport is activated and T peaks below T^n l4l. 
Moreover, at the highest H, their G(T) assumes the form 
for transport by single electrons near the metal to insulator 
transition] 3]. Unlike the NHC films, however, their G(H) 
appear to saturate at 1/Rq = e 2 /h, in the zero T limitl5l 
suggesting a high field quantum metal normal state] 28 ] . This 
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difference could result from stronger spin-orbit scattering in 
Bi NHC films, which is known to quench the quantum metal 
state l29l . Overall, the many resemblances intimate that the 
giant positive MR in the unpatterned films reveals an un- 
derlying CPI phase [Q31 \15\ |T6]| and the subsequent negative 
MR reflects the destruction of pairing by H. Previously, 
the strongest evidence for this association was based on Hall 
measurements! 30] and scaling analyses supporting the dirty 
boson picture fT3ll . 

A notable and perhaps illuminating difference between the 
NHC and the unpatterned films is the extra structure in the 
NHC films' To(H). Inspection of Fig. [3] suggests that 
T = T^ eak (H) + T£ SC (H), where T^ eak (H) is a peaked 
function similar to that exhibited by the unpatterned films and 
Tq sc (H) is a decaying, oscillatory function. The latter likely 
originates through extra terms in the Hamiltonian of the form 
—EjQCos((j)i — <j)j — 2irAij/$o), where the 0's are phases 
on neighboring nodes of the array, is the line integral 
of the vector potential between nodes, and Ejo is the zero 
field phase coupling energy. Such terms are thought to pro- 
duce oscillations of the activation energy of insulating micro- 
fabricated Josephson Junction arrays ifTvTl and can also account 
for decaying oscillations of T c0 of geometrically disordered 
wire arrays (26). The absence of oscillations in the unpat- 
terned films does not necessarily imply the absence of a sim- 
ilar Josephson term in their Hamiltonian. Rather, the absence 
might reflect strong disorder in the geometry of an effective 
wire array model of the films EH EH. The origin of Tfi eak (H) 
is less clear. The disappearance of oscillations near the peak, 
the negative MR and the recovery of clear signatures of sin- 
gle electron transport at high fields could indicate the onset of 
quasiparticle dominated transport at the peak. Consistently, 
the close correspondence between H p and H C 2 suggests that 
TQ Cak (Hys rise at low fields involves the suppression of the 
order parameter amplitude. 

In summary, we have shown that the magnetoresistance of 
amorphous nano-honeycomb Bi films exhibits a giant, posi- 
tive peak, which is similar to that observed in other unpat- 
terned, ultrathin film systems near the superconductor to in- 
sulator transition. The main result is the verification that this 
peak arises from the low field positive magnetoresistance of 
a Cooper Pair Insulator phase with transport dominated by 
incoherent tunneling of Cooper pairs and a high field nega- 
tive magnetoresistance associated with the destruction of the 
pairs. While these conclusions agree qualitatively with some 
models, fundamental questions about the processes driving 
the Cooper pair localization, the emergence of the Cooper Pair 
Insulator phase at high disorder and the size of the localized 
states require further experimental and theoretical attention. 
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